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This  report  presents  two  aspects  of  e  program  designed  to  calculate  elec¬ 
tromagnetic  scattering  from  rough  terrain;  (1)  the  use  of  statistical  estimation 
techniques  to  determine  topographic  parameters  and  (2)  the  results  of  a  single¬ 
roughness -scale  scattering  calculation  based  on  those  parameters,  including 
comparison  with  experimental  data. 

In  the  statistical  part  of  the  present  calculation,  digitized  topographic 
maps  are  used  to  generate  data  bases  tor  the  required  scattering  cells.  The 
application  of  estimation  theory  to  the  data  leads  to  the  specification  of  «--q> 
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i  statistical  parameters  for  each  cell.  The  estimated  parameters  are  then  used 
in  a  hypothesis  test  to  decide  on  a  probability  density  function  (PDF)  that 
represents  the  height  distribution  in  the  cell.  Initially,  the  formulation  uses 
a  single  observation  of  the  multivariate  data.  A  subsequent  approach  involves 
multiple  observations  of  the  heights  on  a  bivariate  basis,  and  further  refine¬ 
ments  are  being  considered. 

The  electromagnetic  scattering  analysis,  the  second  topic,  calculates  the 
amount  of  specular  and  diffuse  multipath  power  reaching  a  monopulse  receiver 
from  a  pulsed  beacon  positioned  over  a  rough  earth.  The  program  allows  for 
spatial  inhomogeneities  and  multiple  specular  reflection  points.  The  analysis 
of  shadowing  by  the  rough  surface  has  been  extended  to  the  case  where  the 
surface  heights  are  distributed  exponentially.  The  calculated  loss  of  bore - 
sight  pointing  accuracy  attributable  to  diffuse  multipath  is  then  compared  with 
the  experimental  results.  The  extent  of  the  specular  region,  the  use  of 
localized  height  variations,  and  the  effect  of  the  azimuthal  variation  in  power 
pattern  are  all  assessed. 
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Electromagnetic  Wave  Scattering  From  Rough  Terrain 


1.  INTRODUCTION 

This  determination  of  the  electromagnetic  scattering  from  rough  terrain  is 
divided  into  two  aspects:  the  statistical  analysis  of  the  terrain  features  and  the 
related  electromagnetic  calculation.  The  statistical  part  involves  specification  of 
terrain  parameters  by  the  use  of  estimation  theory  and  characterization  of  the  topo¬ 
graphic  heights  by  probability  density  functions  (PDFs).  The  electromagnetic  part 
applies  some  of  the  statistical  results  in  a  single -roughness -scale  scattering  cal¬ 
culation  that  will  be  compared  to  experimental  data. 

The  characteristics  of  electromagnetic  signals  scattered  from  rough  terrain 

1  2 

include  contributions  from  clutter  return  and  multipath  return.  ’  These  two 
aspects  can  be  described  by  the  theory  of  scattering  from  rough  surfaces  if  proper¬ 
ties  of  the  terrain  such  as  the  PDF  for  the  surface  height  distribution,  the  covariance 

2 

matrix,  R  the  variance  in  surface  height,  o  and  the  complex  dielectric  constant 

characterizing  the  surface  are  known.  The  numerous  theoretical  models  of  EM 

1  -5 

wave  scattering  from  rough  surfaces  all  relate  the  normalized  cross  section  of 
terrain  to  Iho  foregoing  parameters  characterizing  the  rough  surface. 

In  our  case,  the  physical  parameters  of  the  rough  surface  are  obtained  from 
digitized  terrain  maps  (furnished  by  the  Electromagnetic  Compatibility  Analysis 

(Received  for  publication  20  October  1980) 

f I iii**  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 

So-  References,  page  02.) 
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Center,  ECAC,  and  the  Defense  Mapping  Agency.  DMA).  Estimation  theory  is 

employed  to  specify  the  corresponding  statistical  parameters,  A  hvpothesis  testing 

procedure  determines  the  PDF  for  the  surface  heights. 

The  specific  problem  used  as  an  example  is  that  of  characterising  a  large 

2 

terrain  region  considered  to  be  made  up  of  smaller  aubc “eas  (~4  krr.  ).  The  main 
feature  of  interest  is  the  distribution  of  heights  within  these  subregions.  Each 
subregion  is  characterized  by  a  geologic  code  and  several  statistical  parameters. 

In  particular,  we  are  concerned  with  being  able  to  associate  a  PDF  with  the  range 
of  heights  (z.)  in  the  subregions  and  to  determine  parameters  that  make  the  general 

PDF  explicit.  The  data  elements  z.  *  z^  (x^.  y^l,  where  i  =  1,  2,  3 . N,  and 

N  is  the  total  number  of  grid  points  in  the  x-y  plane  constituting  the  subregion. 

Here,  x^  denotes  the  equally  spared  x-value  along  the  x-axis  and  y  denotes  the 

f**1  y-value  along  the  y-axis,  where  k  =  1,  2 .  and  f  =  1,2 . .  ^*N. 

Thus,  the  N  points  are  distributed  In  the  x-y  plane  so  as  to  form  a  rectangular  grid. 
In  this  analysis,  the  covariance  matrix  can  be  assumed  to  have  the  torm: 

R  =  a2  exp<-T2  /T2]  (1) 

mn  t  mn  ) 

where 

T  =  correlation  Length 


and 


9  2 

T  =  (x  -  X  ) 

mn  m  n 


(v  -  v  r 
m  •  n 


for  the  particular  class  of  data  sets  used. 

Motivation  for  assuming  a  covariance  matrix  of  this  form  is  that  it  leads  to  a 
tractable  mathematical  expression  for  the  incoherent  power  that  is  scattered  when 
an  electromagnetic  wave  is  reflected  from  a  rough  surface.  Once  the  data  have 
been  used  to  specify  the  various  parameters  associated  with  possible  distributions 
of  heights  in  a  region,  several  different  methods  have  been  developed  to  distinguish 
between  distributions  by  means  of  a  binary  decision  hypothesis  testing  procedure. 
In  particular,  the  decisions  involve  whether  certain  Gaussian  or  exponential  dis¬ 
tributions  are  more  appropriate  for  the  observed  heights.  This  specialization 

also  is  motivated  bv  the  theorv  of  electromagnetic  wave  scattering  from  rough 

,1.2 

surfaces. 
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The  normalised  radar  cross  section  of  the  rough  surface  used  in  this  study  was 
6  7  6 

derived  by  Hagfors,  Barrick  and  Semenov.  This  cross  section,  is  incor¬ 
porated  into  a  computer  program  that  calculates  the  amount  of  specular  and  diffuse 
multipath  power  entering  a  monopulse  receiver  from  a  beacon  located  over  rough 
terrain  (see  Figure  1).  The  computer  program  also  calculates  the  error  in  bore- 
sight  pointing  accuracy  of  a  monopulse  receiving  antenna  due  to  noise  and  diffuse 
multipath.  The  computer  program  takes  into  consideration,  among  other  things: 

(1)  the  spatial  nonuniformity  of  the  rough  earth  (that  is,  the  preceding  characteriza¬ 
tion  parameters).  (2)  nonuniformities  in  the  glistening  surface,  (3)  finite  pulse 
length  of  the  beacon,  (4)  antenna  elevation  power  pattern  of  the  monopulse  receiver. 
(5)  multiple  specular  reflection  points  due  to  unevenness  in  surface  heights.  (6)  in¬ 
terference  between  direct  signal  and  multiple,  specularly  reflected  signals,  and 
(7)  finite  azimuthal  beamwidth  of  transmitter  and  receiver.  Finally,  Tor  the  case 

9 

of  normally  distributed  surface  heights.  Sancer's  results  are  employed  to  des¬ 
cribe  the  effects  of  shadowing,  and  for  the  exponentially  distributed  case.  Brown's10 
backscatter  shadowing  calculation  has  been  extended  to  handle  the  forward -scatter 
shadowing  situation. 

The  results  of  the  program  describe  the  effect  of  the  terrain  on  the  electro¬ 
magnetic  signal.  The  final  outputs  include  total  coherent  and  diffuse  power  levels 
and  the  induced  boresight  error.  The  data  output  from  the  computer  program  for 
the  sum  pattern  coherent  power  and  angular  error  in  boresight  is  compared  with 
experimental  data  taken  by  personnel  at  the  MIT  Lincoln  Laboratory  Discrete 
Address  Beacon  System  (DABS)  test  site.  11 


6.  Hagfors,  T.  ( 1964)  Backscattering  from  an  undulating  surface  with  applications 

to  radar  returns  from  the  Moon,  J.  Geophya.  Res.  6^:3779. 

7.  Barrick,  D.  E.  (1968)  Relationship  between  slope  probability  density  function 

and  the  physical  optics  integral  in  rough  surface  scattering,  Proc.  IEEE 
58:1728. 

8.  Semenov,  B,  I.  (1965)  Scattering  of  electromagnetic  waves  from  restricted 

portions  of  rough  surfaces  with  finite  conductivity,  Radiotekh,  i  Elektron 
10:1952. 

MM 

9.  Sanccr,  M.  1.  (1969)  Shadow-corrected  electromagnetic  scattering  from  a 

randomly  rough  surface,  IEEE  Trans,  on  Antennas  and  Prop. 

APjl7;577 -585. 

10.  Brown.  G.S.  (1980)  Shadowing  by  non-Gaussian  random  surfaces.  Proceedings 

of  the  Second  Workshop  on  Terrain  and  Sea  Scatter.  George  Washington 
University,  Washington,  D.  C. 

11.  McGarty,  T.P.  (197  5)  The  Statistical  Characteristics  of  Diffuse  Multipath  and 

its  Effect  on  Antenna  Performance.  AD-A009369. 


9 


Figure  1.  Boxea  Along 
Trajectory  Which  Contribute 
to  Specular  and  Diffuse 
Multipath 


2.  STATISTICAL  ANALYSIS 

For  the  statistical  analysis  of  terrain  heights  in  a  region,  our  approach  has 
been  quite  general.  Depending  on  mission  or  system  requirements,  different 
approaches  to  such  characterization  may  be  more  appropriate.  Two  specific 
formulations  will  be  described;  additional  techniques  are  being  investigated  for 
future  application.  The  two  methods  which  we  have  formulated  in  detail  are  first, 
a  single -observation  multivariate -relation  case  and  second,  a  bivariate  time  series 
approach.  The  results  obtained  by  applying  the  first  approach  to  a  terrain  data 
base  will  be  presented  here.  The  second  is  currently  being  applied  to  the  same  data 
base  and  those  results  eventually  will  be  available  for  comparison. 

2. 1  Single  Obeervation  Formulation 

The  first  approach  involves  a  number  of  operations.  At  the  outset,  we  propose 
some  multivariate  probability  density  functions  that  may  represent  the  observed 
distribution  of  height  values.  Next,  we  use  the  available  data  to  generate  appro¬ 
priate  estimators  of  the  parameters  of  the  respective  densities,  rnally,  we  con¬ 
duct  a  hypothesis  test  to  ascertain  which  density  function  is  more  likely  to  have 
produced  the  observed  height  data. 

The  two  PDF's  for  the  heights  are  multivariate  Gaussian  or  exponential.  The 

12 

Gaussian  has  a  well-known  form: 

12.  Mood,  A.  M. .  and  Graybill,  F.A.  (1963)  Introduction  to  the  Theory  of 
Statistics.  McGraw-Hill.  ~~ 
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- 1 


p(zj.  Zj.  ....  2^  >  =  l(?ir)N^2  I  r|  '  ^2]  exp(--j  (z  -  R  1  (**  -  M 


(2) 


where  R  represents  the  covariance  matrix.  For  our  case,  we  assume  equal  means 

2  2 
(H)  and  equal  variances  (a  )  and  that  the  correlation  function  (p^  =  .Jo  )  has  a 

Gaussian  dependence  on  the  separation  between  points.  The  next  aspect  is  the 

development  of  a  similar  form  for  the  exponential: 


PE(z'. 


z^>  =  Cj  exp(-C2((z'-y)T  R* 1  (z '  -  y  l)1 ' 2) 


■  2, 


(3) 


The  two  coefficients  (C^.  C2)  have  to  be  determined.  This  will  result  in  a  form 
that  satisfies  the  requirements  for  a  PDF.  To  evaluate  the  coefficients,  we  use 
the  properties  that  the  zeroth  moment  integral  of  a  PDF  is  equal  to  unity  and  that 
the  second  moment  integral  is  equal  to  the  variance.  iVe  thus  obtain  a  form  for 
the  multivariate  exponential: 


N*  1 


PE<4 


'2* 


V" 


(2  T.) 


!  1/2-1 


I  ;n  +  i) 


X  exp  (-  ((N*  1)  (z'  -  ix>T  R*1  <z'  -  (x)]1/2  )  .  (4) 


In  order  to  decide  which  of  these  two  PDF's  is  more  appropriate  for  the  data, 
we  must  next  establish  estimators  for  parameters  of  the  densities:  the  means, 
variance,  and  covariances.  The  complexity  of  correlated  multivariate  analysis 
and  the  computational  limitations  associated  with  the  available  data  formats  caused 
us  to  select  estimrtors  that  have  an  intuitive  appeal  based  on  their  form,  rather 
than  a  rigorous  derivation.  From  our  assumption  of  equal  means  and  variances, 
we  use  the  sample  mean  as  the  estimator  for  the  mean  height  and  the  sample 


variance  as  the  variance 

estimate: 

Vn 

Vn 

7  =  (1/N)  E 

E 

z  Uj.yj) 

(mean)  and 

(5) 

1=  1 

j=  i 

Vn 

Vn 

o2  =  (1/N)  E 

E 

<■„ '  ^ 

(variance)  . 

(6) 

i  -  1 

i- 1 

The  procedure  for 

the  covariance  matrix  estimators  is  more  complicated 

.  A 

correlation  length. 

T,  is  defined  as  the  separation  at  which  a  normalized 

covariance 

- 1  2 

function  C  had  decreased  to  the  value  e  ,  where  C  =  y  /u  ,  and  y  is  the 

22  ZZ  ZZ  ZZ 

estimator  for  the  covariance.  The  data  is  used  to  determine  the  estimate  of  T  in 
this  fashion  and  then  the  complete  covariance  estimator  is  formed  from  the  relation 
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-  ff2  exp(-r^n/T2)  (7) 

(separate  x  and  y  relations  are  calculated).  The  form  of  the  covariance  (x  direction) 
used  to  determine  T  is 


Cz2(k)>(l/Nff‘) 


‘  >/n  Vn-k 

I  I 

.  j=  1  i  -  1 


zij  Zi+M 


yfif  >/N-k 

I  I 

3=  1  i«  1 


where  k  represents  the  separation  distance.  A  least-squares-fit  of  the  Czz  values 

to  a  parabola  is  then  used  to  find  T.  Under  the  assumptions  that  have  been  made. 

13 

the  above  estimators  are  similar  to  those  found  in  Jenkins  and  Watts.  Values 

for  T  in  the  x  and  y  directions  were  obtained,  and  then  the  two  values  averaged  to 

find  the  final  estimate  of  T  in  a  given  subregion. 

We  now  have  the  two  PDF  forms  and  the  required  parameters.  The  next 

aspect  is  to  decide  which  PDF  is  the  more  appropriate  for  the  given  data.  One 

final  comment  on  the  parameters:  In  order  to  satisfy  the  restrictions  of  the  PDF's, 

it  is  necessary  to  show  that  the  quadratic  form  appearing  in  both  cases  is  positive 

definite.  This  has  been  demonstrated  for  the  foregoing  cases  by  making  use  of 

the  Gaussian  form  assigned  to  the  covariance  matrix  elements.  Details  can  be 

1 4 

found  in  the  report  by  Lennon  and  Papa. 

The  form  of  the  hypothesis  test  used  here  is  based  on  the  maximum  a  posteriori 
probability  criterion.  This  is  equivalent  to  a  minimum  error  probability  criterion. 
We  assign  hypothesis  Hj  to  the  Gaussian  case  and  hypothesis  Hq  to  the  exponential. 
Then  the  likelihood  ratio  parameter. 


A  =  Pj  (z j.  a 2 


ZN*  !  Po  *Z1'  Z9 . ZM*  • 


Let  P(Hq)  be  the  probability  that  hypothesis  Hq  is  true,  Then  the  decision  rule  may 
be  written  as:  Choose  Hj  if 

P(H  > 

\  >  a  l 


For  our  case,  we  assume  that  it  is  equally  probable  that  hypothesis  Hj  or  Hq  is 
true  and  the  decision  rule  reduces  to  whether  or  not  A  -  1.  Note  that  it  may  be 


13.  Jenkins,  G.M.,  and  Watts,  D.G.  (  1968)  Spectral  Analysis  and  Its  Applications, 

Holden -Day,  * 

14.  Lennon,  J,  F. ,  and  Papa,  R.J.  (1980)  Statistical  Characterization  of  Rough 

Terrain,  RADC-TR-80-9,  RADC/EE  Hanscom  AFB,  Massachusetts. 
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possible  to  alter  the  probability  that  Hq  is  true  based  on  external  evidence  (such 
as  the  type  of  terrain). 

When  the  specific  forms  for  the  two  PDP's  are  introduced  into  this  relation 
it  becomes 


(11) 


where 


C*  =  |(z  -  y)T  R*1  (z  -  fct)] 1/2  . 


For  convenience,  we  rewrite  the  test  in  logarithmic  form  and  obtain  the  result  that 

+  (N  +  1>  -  In  3  -  N  In  1/2 


Hj  is  true  if  either 


Q  S 


-(rm) 


+  \IN+  1 


(12) 


or 


Q  ^  j  yfsTl  -  ^2  ln(r  (N+  1)  -  In  ir  -  N  ln(^iyl)|1/2  .  (13) 


For  the  actual  cases,  N  «  100,  and  the  specific  result  is  our  decision  that  the 
terrain  heights  in  a  given  region  are  from  a  Gaussian  PDF  if 

85.01  S  [<z  -  £>T  R*1  <z  -  ji)]  S  119.37 

2  2 

and  conversely  the  points  are  exponential  if  Q  >  118.  37  or  Q  <  85.  Cl. 


2.2  Multiple  Ob*erv*tion  Formulation 

The  second  approach  to  statistical  characterization  of  terrain  heights  eliminates 
the  need  for  matrix  inversion  by  dealing  with  the  points  as  successive  sets  of  pairs. 
This  allows  us  to  use  a  bivariate  form  of  the  probability  density  function  (PDF) 
when  constructing  the  hypothesis  test.  This  new  approach  is  aimed  at  characteriz¬ 
ing  the  relations  between  pairs  of  points  in  a  region  where  we  assume  isotropy  and 
atatlonarlty  In  the  height  distribution.  This  Is  an  appropriate  characterization, 
since  the  related  electromagnetic  calculation  Is  concerned  with  the  mean  value  of 
the  scattered  power. 

The  two-point  characterization  is  a  supplement  to  the  original  N-variate  result; 
it  is  from  the  standpoint  of  calculation  more  appealing  and  should  give  better  results 
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by  avoiding  the  problems  of  inverting  large  matrices  and  by  introducing  a  form 
that  allows  multiple  observations,  rather  than  a  single  multivariate  observation. 

As  in  the  foregoing  discussion  two  basic  types  of  analysis  are  involved:  (1) 
he  determination  of  appropriate  estimators;  and  (2)  a  hypothesis  test  using  PDF's 
based  on  the  estimators.  This  introduces  one  complication;  that  is,  the  estimators 
are  always  based  on  summations  Involving  N  =  10,  whereas  the  test  uses  subsets 
of  points  in  sums  with  two  limits,  K  ■  48  and  S  =  16;  K  is  the  number  of  different 
pair  spacings  to  be  considered  and  S  the  number  of  observations  at  each  spacing. 

The  estimators  for  the  means  and  variances  are  the  same  as  those  described 
in  the  first  approach.  For  the  covariance  estimation,  however,  the  situation  is 
different.  We  now  form  a  distinct  estimator  for  each  separation:  C  (r),  where  r 
is  an  index  assigned  to  the  combined  (a,b)  grid  spacing. 

Numerical  values  for  Indices  in  the  following  sections  correspond  to  the  specific 

case  where  the  method  has  been  applied  to  the  ECAC  data  base.  The  values  of  r 

are  r  *  0,  1,  2,  ....  K  =  48;  C<r)  is  assigned  to  C(a,  b>  according  to  the  scheme 

described  below,  but  this  is  arbitrary  and  other  orderings  could  be  used.  There 

are  thus  49  covariance  estimators  since  the  x  and  y  coordinate  spacings  on  the  grid 

differ  and  we  allow  a  =  0,  1,  2,  3,  4,  5,  6  and  b>  0,  1,  2,  3,  4,  5,  6  where  a  A  x 

Is  the  x  coordinate  separation  and  b  A  y  is  the  y  coordinate  separation  for  the  two 

points  of  the  pair.  The  number  of  possible  pairs  at  a  given  separation  (a,  b)  is 

(N-a)  x  (N-b).  Thus  for  a  *  0.  b  =  0,  there  are  100  possible  pairs,  (z  ,  z  ) 

pq  p+a,  q+b 

and  for  a  =  6,  b  =  6,  there  are  16  possibilities.  Then  by  assigning  each  of  the 
49  r-values  to  a  particular  (a,  b)  combination  we  have: 


C(r)  =  C(a.b)  = 


.  / (N*a)(N-b)  z2 
+  - 


N-a 

N-b  \ 

L 

p=  i 

YL  z  z  .  ) 

q=i  pq  FH-a,  q+b J 

I  N-a  N-b  v* 

t 

/  £  (z  +  z  .)  | 

l  p.  1  qM  ™  *»•'  q+b  j 

(14) 

..2  2 
N  <7 

where  for  convenience  we  assign: 


C(0)  •*  C(0, 0); 

C(l) 

"  C(0.  1); 

C<6)  - 

C(0,  6); 

C(7)  **  C(l,  0); 

C(8) 

C(27) 

-  CU,  1); 

-  C<3,6); 

#  •  • 

C(13)  "* 

C(l,  6); 

C(4i)  **  C(5,  6); 

and 

C(48)  -* 

C(6,  6). 

(Note  that  with  this  notation  C(0)  =  C(0,  0)  =  1.  This  is  consistent. ) 
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At  this  point  we  have  al]  the  necessary  estimators  for  use  in  the  hypothesis 
test  and  we  could  proceed  directly  to  that  aspect.  However,  for  the  electromag¬ 
netics,  it  U  still  desirable  to  determine  the  correlation  length  T,  in  a  similar 

fashion  to  the  original  case,  although  it  will  not  be  used  in  the  statistics.  We 

2  2  2 

already  have  our  C(r)  values  defined,  so  we  set  d  (r)  *  (a  Ax)  +  (b  A  y)  and 
obtain  as  before: 


48 

£ 

r=  1 


d2<r)  - 


d4(r) 

TB - 

£  d2(r)  C(r) 

r  =  1  ' 
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in  meters. 


(15) 


The  nt  1  topic  is  that  of  the  hypothesis  test.  Here,  as  we  pointed  out  before, 
each  r  value  will  use  only  S  =  16  cases.  For  (a  »  6,  b  =  6)  there  are  only  16 
possible  pairs  of  (s  ,  Zp+a  q+b^’  but  *or  the  other  cases  there  are  additional 
possibilities.  Some  scheme  has  to  be  used  [that  is,  random  number  generation,  or 
selection  from  a  broad  range  of  (p,  q)  values]  so  that  a  representative  selection  of 
pairs  are  used  in  the  test  for  each  r  value. 

The  hypothesis  test  is  similar  to  the  previous  form  but  in  this  case  we  have 
multiple  observations  and  so  the  probabilities  with  which  we  are  concerned  are  the 
likelihood  functions  for  the  respective  bivariate  probability  densities.  Since  we  are 
dealing  with  multiple  observations,  it  should  be  noted  that  in  the  following  relations, 
the  superscripts  identify  the  two  zp  ^  or  z  a  q+b  members  of  the  pair,  whereas 
the  subscripts  r  and  i  now  refer  to  the  separation  and  selection  number,  respectively. 

For  the  Gaussian  case 

p,  i-y  jKS  J(i  -  c2a^(i  -  c2<2>)  (i  -  c2(k»)J  7 


(16) 


For  the  exponential  case 


s 


i  KS  f 

p°‘(wr)  1 

C2<r))  (l 

■  C2!2>)  • 

(l  -  C*(K>)j  2 

B 

X  s 

r  t 

aon  -i)M?  -).<•!? 

*  exp 

r-  1  1-  l' 

^  1  -  c2(r) 

(17) 


15 


These  terms  are  then  introduced  into  the  hypothesis  test  and  simplified  by¬ 
formation  of  the  log -likelihood  function,  which  leads  to  the  following  decision: 
Choose  a  Gaussian  PDF  when  TEST  ~  S  K  (3  -2  In  3)  a  ^ ,  This  becomes 


i 


48 

£  G(r) 

r  =  l 


-  6170  S  0 


where 


G(r) 


16 

E 

i=  i 


f/_U)  -2  ,J1)  ru,(2)  r\. /_(2) 

'  ri  "Z)  -2C(r)(*rl  -*)(*rl  "z)+(zri  '  z“ 


1  -  C  (r) 


l/2  2 

)  ->F 

(18) 


The  results  for  this  approach  are  not  yet  available.  The  single  observation, 
r  lultivariate  characterization  is  thus  the  only  one  that  can  be  used  as  an  input  to 
the  electromagnetic  calculations. 


3.  SITE  CHARACTERIZATION 

The  specific  site  used  for  the  characterization  is  one  in  eastern  Massachusetts.  11 
This  area  was  selected  because  of  the  availability  of  electromagnetic  terrain 
scattering  data  obtained  during  tests  of  a  beacon  system  at  that  location.  A  rec  - 
tangular  area  around  this  Discrete  Address  Beacon  System  (DABS)  site  was  desig¬ 
nated.  The  rectangular  area  was  43.3  km  long  and  42.8  km  wide.  The  area  was 
then  subdivided  Into  smaller  rectangular  cells,  each  with  sides  of  2050  m  by 
1825  m.  Each  cell  is  further  subdivided  into  a  10  by  10  grid  of  points.  A  data  base 
of  topographic  elevations  for  this  area  is  available  at  the  Electromagnetic  Compati¬ 
bility  Analysis  Center  (ECAC).  This  was  prepared  from  Defense  Mapping  Agency 
(DMA)  supplied  digitized  terrain  maps  at.  1:  250,000  scale  size. 

The  statistical  data  for  each  ceil  has  been  recorded  on  a  computer  tape  for  use 
with  the  program  for  the  electromagnetic  analysis.  Each  cell  is  represented  by 
seven  descriptors.  The  first  two  entries  are  the  <x,y)  coordinates  for  the  center 
of  the  cell.  (The  origin  of  the  coordinate  system  is  taken  as  the  center  of  the  ex¬ 
treme  southwestern  corner  of  the  rectangular  region.  )  The  next  item  is  the  geo¬ 
logical  code  for  the  cell.  The  predominant  feature  is  woods;  there  are  a  nun. her  of 
cells  containing  clusters  of  lakes  and  ponds  and  a  few  town  sites  with  associated 
cleared  areas.  This  is  followed  by  the  mean  and  variance  of  the  heights  in  the  cell 
and  the  estimated  correlation  length,  Y  (the  units  of  length  are  in  meters).  The 
final  quantity  is  the  result  of  the  hypothesis  test.  This  result  is  presented  in  a 
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format  such  that  the  heights  in  a  region  are  Gaussian  when  (-S3.  36  S  TEST  5  0). 
When  TEST  *  0  and  the  variance  is  very  small,  the  region  is  essentially  a  smooth 
surface  (no  roughness). 

14 

When  the  over-all  results  for  the  region  are  examined,  one  observation  that 

can  be  made  is  that  when  the  magnitude  of  TEST  ta  very  large,  the  correlation 

length  is  also  very  large;  for  those  cases  T  is  comparable  to  one  half  the  cell  size 

or  even  larger.  When  that  occurs,  the  determinant  of  the  covariance  matrix,  R, 

becomes  very  small.  As  a  result.  It  is  increasingly  difficult  to  obtain  an  accurate 

Inverse  of  g,  due  to  rapid  build-up  of  round-off  error.  It  should  be  noted  here  that 

the  second  characterization  technique  is  based  on  bivariate  formulation,  so  the 

Inverse  of  R  is  introduced  analytically.  Hence,  in  that  case,  the  round-off  errors 

of  the  machine-calculated  inverse  do  not  enter  the  result.  A  related  difficulty  in 

the  results  of  the  first  approach  can  be  seen  in  those  cases  where  TEST  <  -  118.  37. 

2 

This  would  be  possible  only  if  the  quadratic  form  Q  is  not  positive  definite.  This 
contradiction  of  that  theoretically  imposed  requirement  implies  that  further 
machine-induced  errors  were  present  and  those  results  can  not  be  considered  valid. 

The  grid  structure  for  the  site  in  Massachusetts  has  been  analyzed  in  terms  of 
specific  trajectory  of  the  test  program  in  an  attempt  to  obtain  a  quick  correction 
of  the  hypothesis  test  results.  The  actual  subset  of  boxes  that  contribute  to  the 
specular  and  diffuse  scattering  analysis  for  the  trajectory  are  shown  in  Figure  1. 
For  those  cases,  the  analysis  was  reproduced  on  a  more  accurate  computer  and  a 
second  set  of  TEST  results  were  obtained.  This  second  set  did  not  suffer  from  the 
extreme  round-off  problems  of  the  original  calculation.  These  results  have  been 
used  in  the  present  comparisons. 

In  order  to  use  these  results  in  the  rough  surface  electromagnetic  calculations, 

one  additional  aspect  should  be  noted.  For  the  types  of  geological  features  that 

describe  the  respective  regions,  data  exist  on  the  associated  complex  dielectric 

15  3 

constants  at  microwave  frequencies.  ' 

4.  ELECTROMAGNETIC  SCATTERING  ANALYSIS 

The  radar  cross  section  of  terrain  is  normalized  with  respect  to  the  average 
area  illuminated  by  the  radar.  The  normalized  cross  section,  0 Q.  may  be  divided 


14.  Lennon,  J.  F. ,  and  Papa,  R.  J.  (1980)  Statistical  Characterization  of  Rough 

Terrain,  RADC-TR-80-9,  RADC/EE  Hanscom  aPb,  Massachusetts. 

15.  Lytle,  R.  J.  { 1 9 <  4 )  Measurement  of  earth  medium  electrical  characteristics? 

Techniques,  Results  and  Applications,  IEEE  Trans,  on  Geoscience 

Electronics,  GE-12-.81  “ 

■  ■  ■  1 


17 


16 

into  three  general  categories:  (1)  the  slightly  rough  surface;  (2)  the  very  rough 
surface,  7'  ®  and  (3)  the  multiple  scale  rough  surface. This  study  deals 
with  the  second  category,  that  is,  the  very  rough  surface  (irregularities  are  large 
compared  to  a  wavelength).  Fine  scale  effects  are  introduced  cursorily  by  the  use 
of  effective  complex  dielectric  constants  assigned  on  the  basis  of  geologic  des¬ 
criptors  in  each  terrain  region  (lakes,  woods,  roads,  and  so  or.).  This  is  an 
adequate  model  for  the  forward  scattering  problem  associated  with  beacon  trans¬ 
mission,  but  for  cases  where  backscatter  results  are  of  interest  the  fine  scale 
contributions  have  to  be  considered  in  greater  detail. 

This  is  justified  for  the  following  reason.  At  low  incident  grazing  angles,  the 
large,  gently  undulating  irregularities  will  tend  to  scatter  in  the  forward  direction 
and  represent  the  dominant  contribution  in  that  direction.  The  fine  scale  rough¬ 
ness  (where  irregularities  are  small  compared  to  a  wavelength)  will  tend  to  scatter 
energy  in  directions  satisfying  the  Bragg  conditions,  as  modified  by  the  large  scale 
tilting  effects.  These  fine  scale  results  will  dominate  in  the  backscatter  region, 
where  there  are  only  secondary  large  scale  roughness  contributions,  but  only 
represent  a  secondary  contribution  to  the  complete  forward  scatter  result. 


4.1  Scattering  Ctom  Section 

2 

Ruck  et  al  give  expressions  for  the  average  bistatic  rough  surface  cross 

section  a  under  the  following  four  assumptions:  (1)  the  radius  of  curvature  of  the 

surface  irregularities  is  larger  than  a  wavelength;  (2)  the  roughness  is  isotropic 

in  both  surface  dimensions;  (3)  the  correlation  length  is  smaller  than  either  the 

x  or  y  dimension  of  the  sample  subregion;  and  (4)  multiple  scattering  is  neglected 

2 

Using  the  notation  of  Ruck  et  al,  one  finds  that  the  expression  for  becomes 

*o=l/3pq|2J  (19) 


where 

J  =  (T2/o24  2  )  exp 


+  « 

T~ 
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for  a  Gaussian  bivariate  surface  height  probability  density  function  and 


J=  (3T 2lo2l2z)  exp 


M(¥T1 


(20) 


for  an  exponential  surface  height  probability  density  function.  The  scattering 
matrix  elements  are  given  by 


hh 
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£  =  -cos  9 .  ~  1.0 
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f).  =  ang'e  of  incidence  (with  respect  to  surface  normal) 


=  angle  of  scattering  (with  respect  to  surface  normal) 


and 


Here,  e  is  the  relative  complex  dielectric  constant  of  the  surface,  the  subscript 
||  refers  to  the  E-iield  in  the  plane  of  incidence,  and  the  subscript^  refers  to  the 
E-field  normal  to  the  plane  of  incidence.  These  simolified  forms  of  Ruck's  ex¬ 
pressions  follow  from  the  assumption  that  the  receiver  is  far  from  the  transmitter 
so  that  the  portion  of  the  "glistening  surface"  ^  that  contributes  to  the  diffuse  multi - 
path  is  a  long,  narrow  strip  extending  between  transmitter  and  receiver.  This 
assumption  allows  iis  to  make  the  approximation  that  the  azimuthal  scattering  angle, 
9^  ~  0.0.  This  will  be  further  discussed  in  Section  4,5. 


4.2  Shadowing 

The  effect  of  shadowing  on  the  diffuse  scatter  from  the  surface  is  introduced 
by  multiplying  the  expression  for  <7^  by  an  appropriate  shadowing  function,  depend¬ 
ing  on  whether  the  surface  heights  of  the  particular  contributing  subregion  are 

normally  or  exponentially  distributed. 

g 

In  this  study,  Sancer's  results  are  used  to  describe  the  effects  of  shadowing 
when  the  surface  heights  are  normally  distributed.  For  the  situation  where  the 
surface  heights  are  exponentially  distributed.  Brown's^  expressions  for  shadowing 
in  the  case  of  backscattering  have  been  extended  to  include  shadowing  in  the  forward 
scattering  direction.  The  details  of  the  exponential  shadowing  formulation  are  pre¬ 
sented  in  Appendix  A. 

4.3  Program  Initiation 

The  computer  program  is  designed  to  incorporate  the  expression  for  aQ  into 
an  integral  over  the  glistening  surface.  The  program  has  the  capability  of  calcu¬ 
lating  the  coherent  power  (specular  plus  direct)  and  diffuse  multipath  power  reaching 
a  monopulse  receiver  from  a  beacon  over  rough  terrain  (Figure  2).  The  computer 
program  uses  the  previously  described  data  tape  of  the  statistical  parameters  for 
a  particular  site  as  an  input.  Other  input  variables  characterize  the  transmitter, 
the  environmental  aspects,  and  the  receiver.  Transmitter  values  include;  the 
gain  of  the  transmitter,  the  polarization  of  the  transmitted  wave,  the  peak  power 
of  the  transmitted  pulse,  the  pulse  length,  and  the  wavelength  of  the  signal.  Ex¬ 
ternal  inputs  include;  the  complex  dielectric  constant  of  each  type  of  geological 
region,  the  coordinates  (latitude  and  longitude)  of  the  monopulse  receiver,  the 
initial  and  final  position  of  the  aircraft  containing  the  transmitter,  a  parameter  to 
control  the  effects  of  shadowing,  and  the  velocity  of  the  aircraft.  Receiver  data 
include;  the  height  of  the  receiver,  the  bandwidth  of  the  receiver,  the  front-end 
receiver-noise  figure,  the  antenna  gain  for  the  sum  and  difference  patterns  of  the 
monopitlse  receiver,  its  azimuthal  beam  width,  the  sampling  frequency  of  the  re¬ 
ceiver,  the  transmission  line  loss  factor  of  the  cables  connecting  the  antenna  to 
receiver,  and  the  difference  pattern  slope  near  the  boresight  axis. 

From  a  knowledge  of  the  Initial  and  final  positions  of  the  aircraft  and  the  air¬ 
craft  speed,  the  computer  program  first  calculates  the  trajectory,  as  a  function  of 
time.  Then,  from  a  modified  form  of  the  radar  range  equation,  the  electric  field 
intensity  of  the  direct  signal  at  the  receiver  is  calculated  at  fixed  time  intervals 
(sampling  time  is  an  input  variable  to  the  program). 
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PULSED  RADAR  BEACON 


MONOPULSE  RECEIVER 


ROUGH  TERRAIN 

Figure  2.  Reflection  of  Radar  Waves  From  Rough  Terrain 


4.4  Specular  Multipath 

For  the  specular  multipath  rays,  all  possible  specular  reflection  points  be¬ 
tween  the  transmitter  and  receiver  are  determined  from  each  position  of  the  trans¬ 
mitter.  Multiple  specular  reflection  points  due  to  unevenness  in  surface  height 
are  taken  into  consideration.  At  each  possible  point,  the  path  to  both  antenna  and 
transmitter  is  examined  for  ray  blocking  by  the  surface  heights  along  the  trajectory. 
The  appropriate  finite  dielectric  constant  of  the  earth  at  each  specular  point,  the 
antenna  elevation  pattern  (receive)  and  the  surface  roughness  are  also  accounted 
for  in  calculating  the  phase  and  amplitude  of  each  specular  multipath  ray.  At  each 
point  on  the  transmitter's  trajectory,  the  total  coherent  power  for  the  sum  and 
difference  channels  in  the  monopulse  receiver  is  calculated. 

One  aspect  of  the  discussion  of  specular  multipath  that  relates  to  the  magnitude 
of  the  effect  as  a  function  of  range  is  the  extent  of  the  contributing  region.  Figure  3 
shows  the  history  of  the  region  as  the  target  approaches  the  receiving  antenna.  For 
this  analysis,  the  first  Fresnel  zone  is  considered  the  area  from  which  specular 
reflections  can  occur  and  the  remaining  zones  are  assumed  to  cancel  each  other. 

The  curves  depict  the  specular  point  location,  the  center  of  the  Fresnel  zone,  and 
its  extent  and  width.  Details  of  these  calculations  are  included  in  Appendix  B  and 
the  significance  of  the  results  is  analyzed  in  Section  5. 


4.5  Hw  Glutening  Surface 

In  the  calculation  of  the  total  diffuse  power,  account  is  taken  of  spatial  in¬ 
homogeneities  of  the  rough  earth  and  nonuniform ities  in  the  boundaries  of  the 
glistening  surface.  In  Figure  4,  a  typical  illustration  is  given  of  this  surface.  The 
distance  from  the  receiver  to  the  edge  of  the  glistening  surface  is  lj,  and  the  dis¬ 
tance  from  the  transmitter  to  the  opposite  edge  of  the  surface  is  denoted  fg*  These 
distances  are  a  function  of  the  ratio  (7/T,  and  explicit  expressions  for  them  are  given 
by  Beckmann  and  Spizzlchino,  *  For  the  system  and  environmental  parameters  being 
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investigated  in  this  report,  it  was  fpund  from  computational  experiments  that  sig¬ 
nificant  additional  contributions  to  the  diffuse  power  originated  from  regions  beyond 
the  classical  definition  of  the  glistening  surface.  To  allow  for  this,  the  expressions 

for  f  and  were  modified  so  as  to  extend  the  length  of  this  glistening  surface.  In 
1  * 

calculating  the  total  diffuse  power,  the  normalised  cross  section  Oq  is  integrated 
over  the  glistening  surface.  Since  the  terrain  is  inhomogeneous,  the  integration  iB 
accomplished  by  dividing  the  surface  into  small  strips,  as  shown  in  Figure  4.  The 
total  diffuse  power  is  calculated  by  summing  the  contributions  from  each  strip. 


Hi  I-  H«,h- 


Figure  4.  The  Glistening 
Surface 


4.6  Diffiue  Multipath 

The  diffuse  multipath  power,  Ppjpp,  entering  the  receiver  is  obtained  from 
the  equation: 


where 
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■  transmitted  power, 

«  wavelength, 

«  gain  (power)  of  transmitter  in  azimuth, 

*  gain  of  receiver  in  azimuth, 

=  gain  of  transmitter  in  elevation, 

=  gain  of  receiver  in  elevation, 

=  angle  between  boresight  and  point  on  glistening 
surface  for  transmitter, 

■  angle  between  boresight  and  point  on  glistening 
surface  for  receiver, 

-  range  between  transmitter  and  point  on  glistening  surface, 

=  range  between  receiver  and  point  on  glistening  surface, 

=  element  of  area  of  glistening  surface  which  is  illuminated 
by  beacon 

In  the  numerical  integration  of  the  equation  for  Ppjpp*  the  transmitter  is 
assumed  to  have  a  uniform  azimuthal  power  pattern  and  the  azimuthal  variation  in 
<Tq  has  been  assumed  negligible.  Two  different  cases  were  considered  for  the 
azimuthal  power  pattern  variation  of  the  receiving  antenna.  First,  it  was  assumed 
to  be  constant  at  5  dB  below  the  sum  pattern  gain.  This  gain  corresponds  to  3  dB 
below  the  peak  In  the  difference  gain.  In  a  more  refined  analysis,  the  power  pattern 
was  described  by  a  parabolic  dependence  on  azimuthal  angle  and  treated  as  a 
variable  in  the  integration  over  the  power  pattern.  The  details  of  this  azimuthal 
power  pattern  variation  are  presented  in  Appendix  C. 

4.7  Boresight  Error 

To  calculate  the  error  in  boresight  of  the  monopulse  receiver  due  to  diffuse 
multipath  and  receiver  noise,  it  is  assumed  that  the  diffuse  multipath  is  decorrela- 
ted  from  pulse  to  pulse  and  that  the  spectral  width  of  the  diffuse  multipath  power  is 
narrow  compared  to  the  bandwidth  of  the  receiver/ processor.  These  assumptions 
appear  to  be  justified  on  the  basi9  of  rough,  order-of-magnitude  estimates  of  the 
appropriate  parameters.  Under  these  assumptions,  the  total  amount  of  noiselike 
interference  is  given  by 
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NI  5  PDIFF+No'  where 

Nj  =  noiselike  interference  power 

Nq  3  noise  power  from  environment  plus  receiver. 

The  error,  CA>  in  azimuthal  boresight  pointing  accuracy  is  given  by  the 
“  19 

expression  irom  Barton  and  Ward  as  follows: 


where 

6  =  azimuthal  beamwidth, 

D 

STIR  =  ^coh^I  =  signal  to  interference  ratio  in  the 
difference  channel, 

F*coh  =  coherent  power, 

km  =  normalised  monopulse  slope  (obtainable  from  sum 
and  difference  patterns). 

The  output  of  the  computer  program  consists  of  azimuthal  angular  error  in 
boresight  due  to  noise,  the  error  due  to  noise  plus  diffuse  multipath,  the  total  sum 
pattern  coherent  power,  the  total  diffuse  power,  signal -to-noise  ratio,  signal-to- 
interference  ratio,  and  range  from  transmitter  to  receiver.  In  the  conclusions, 
the  computer  output  for  the  analysis  of  a  particular  site  is  compared  with  some 
experimental  data. 

5.  RESULTS  AND  CONCLUSIONS 

Two  theoretical  results,  the  coherent  power  in  the  sum  channel  of  a  monopulse 
antenna.  pcoh*  and  the  standard  deviation  In  boresight  pointing  accuracy,  <Jq,  are 
compared  with  the  experimental  data  of  McGartv,  11  taken  at  the  DABS  site. 

5.1  Experimental  Condition* 

The  DABS  receiver  was  a.'  L-band  rotatable  array.  The  beacon  was  located 
on  a  U-10  aircraft  which  flew  a  number  of  radial  trajectories  toward  and  away  from 
the  monopulse  receiver.  Data  were  recorded  for  about  100  flights  for  different 

19,  Barton,  D.  K.  ,  and  Ward,  H.  R.  ( 1969)  Handbook  of  Radar  Measurement, 
Prentice -Hall, 
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aircraft  heights,  different  radial  flight  trajectories,  and  different  receiver -antenna 
tilt  angles,  ihe  conditions  under  which  the  data  were  taken  are  listed  in  Table  l. 


Table  1.  Experimental  Conditions  for  DABS  Tests 


Front  end  receiver  noise  figure 

3  dB 

Gain  of  monopulse  receiver  (sum  pattern) 

22.  5  dB 

Gain  of  transmitter 

4  dB 

Height  of  receiver 

101  m 

Height  of  transmitter 

1220  m 

Signal  polarization 

vertical 

Peak  transmitter  power 

350  W 

Pulse  length 

20  fi sec 

Azimuthal  beamwidth  (receiver) 

3* 

Wavelength 

0. 275  m 

Transmission  line  loss  factor 

3  dB 

In  order  to  make  the  comparisons,  various  manipulations  were  required.  In 
the  equation  for  Oq,  the  normalized  slope  in  the  difference  pattern  near  boresight, 
k^.  was  obtained  from  graphs  of  the  sum  and  difference  patterns  of  the  monopulse 
antenna.  Within  the  accuracy  of  these  graphs,  it  appears  that  (1.5  -  km  -  1.7), 
Figures  5  and  6  present  experimental  data  with  which  various  theoretical  results 
are  to  be  compared.  Figure  5  is  a  plot  of  the  sum  signal  ipcoj^  versus  range  of 
transmitter  to  receiver  in  nmi.  Figure  6  is  a  plot  of  azimuth  error  (Oq)  versus 
'ange  in  nmi. 

5.2  Coherent  Power  Results 

The  magnitude  of  the  coherent  power  multipath  contribution  from  specularly 
reflected  rays  depends  on  the  relationship  of  the  specular  region  to  the  variance 
and  correlation  length,  as  well  as  the  over-all  size  of  the  specular  region  used  in 
the  calculation.  The  effect  of  these  aspects  on  the  nature  of  the  coherent  results 
is  significant. 

It  should  b'  '-eealled  that  Figure  3  shows  the  extent  of  the  first  Fresnel  zone 
decreasing  as  range  from  transmitter  io  receiver  decreases.  Figure  1  illus¬ 
trates  all  the  boxes  along  the  trajectory  that  contribute  to  the  specular  and  diffuse 
multipath  powei ;  me  actual  boxes  that  contribute  to  the  result  in  a  given  case  de¬ 
pend  on  the  position  of  the  transmitter  along  the  trajectory.  In  particular,  for  the 
range  of  sv3tem  ana  environmental  parameters  under  investigation  in  this  report, 
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AZIMUTH  ERROR  (DEG)  -SUM  SIGNAL  (dBm) 


only  the  three  boxes  closest  to  the  receiving  antenna  (in  the  southeastern  corner  of 
the  total  region)  contribute  to  the  specular  reflections.  Also,  for  each  position 
of  the  transmitter,  there  is  only  one  specular  reflection  point. 


Initially,  the  variance  associated  with  the  distribution  of  heights  over  an  entire 
box  (see  Figure  1)  was  used  in  the  Rayleigh  attenuation  factor  of  the  specular 
multipath  ray.  The  coherent  power  from  the  specularly  reflected  rays  is  small 
in  that  case  and  the  result  is  just  the  direct-ray  (1/R  )  fall  off  of  power  with 
range  R.  This  result  is  not  inconsistent  with  what  would  be  expected  when  the  size 
of  the  actual  contributing  region  (on  the  order  of  the  first  Fresnel  zone)  is  con¬ 
sidered,  along  with  the  associated  assumption  of  small  correlation  length  contained 
in  the  scattering  formulation.  When  the  first  Fresnel  zone's  dimensions  are  con¬ 
siderably  smaller  than  the  box  size  (2  km  by  2  km),  smaller  sample  regions  should 
be  used  in  calculating  the  variances  in  surface  height. 

The  calculation  of  more  localized  variances  associated  with  the  Fresnel  zone 
dimensions  was  then  performed.  This  resulted  In  smaller  values  for  the  variances 
and  increased  specularly  reflected  multipath  power.  Figure  7  shows  the  calculated 
coherent  power  for  the  sum  signal,  together  with  the  actual  data.  The  specular 
contributions  were  calculated  using  variances  based  on  trajectory -centered  sub- 
regions  of  the  three  boxes  closest  to  the  antenna.  As  the  specular  point  shifts 
across  the  three  regions,  the  effect  of  the  different  values  for  the  three  local 
variances  can  be  seen  In  successive  changes  that  appear  in  the  behavior  of  the  sig¬ 
nal.  The  theoretical  results  tend  to  be  about  6  dB  higher  on  the  average;  this  is 
most  likely  due  to  the  fact  that  the  receiver/processor  losses  are  unknown  and 
hence  have  been  neglected.  When  typical  processing  losses  of  -5  dB  are  included 
In  the  analysis,  the  calculated  coherent  power  for  the  sum  signal,  together  with 
the  actual  data,  are  In  remarkable  agreement,  as  may  be  seen  from  Figure  8. 


Figure  7.  Theoretical  Calculations  and 
Experimental  Data;  Sum  Signal  vs  Range 
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Figure  8.  Theoretical  Calculations  and 
Experimental  Data:  Sum  Signal  vs  Range 
(system  losses  1  -5  dB) 


It  should  be  noted  that  the  usual  expression  for  the  Rayleigh  attenuation  factor 
of  the  specularly  reflected  ray  (as  given  by  Beckmann  and  Spitzichino)  was  de¬ 
rived  by  assuming  that  the  correlation  length  T  is  much  smaller  than  the  dimen¬ 
sions  of  the  Fresnel  zone.  For  close  ranges,  the  parameters  investigated  In  this 
report  are  such  that  T  becomes  comparable  to  the  dimensions  of  the  Fresnel  zone. 

In  this  case,  additional  corrections  have  to  be  made  for  the  specularly  reflected 
ray.  This  problem  area  is  currently  being  investigated. 

5.3  Boredgfit  Error  Remit « 

Before  specific  results  are  addressed,  some  general  comments  should  be  made. 
First,  for  the  range  of  parameters  used  in  this  analysis,  the  inclusion  of  shadowing 
in  the  boresight  error  calculation  does  not  introduce  a  significant  change  from  the 
cases  where  shadowing  is  neglected.  Second,  unlike  the  coherent  power  results, 
present  results  extend  to  a  maximum  of  only  30  nmi  (the  trajectory  originates  at 
30  nmi).  The  coherent  results  are  affected  only  by  the  area  close  to  the  receiving 
antenna  whereas  the  glistening  surface  extends  over  a  considerable  extent  of  the 
trajectory  for  the  diffuse  multipath  power.  Therefore,  the  finite  nature  of  the 
geographical  data  base  limited  the  farthest  point  at  which  the  diffuse  contributions 
could  be  assessed.  Finally,  it  should  be  noted  that  these  results  do  not  include 
any  fine  scale  roughness  contributions  to  the  diffuse  power. 

The  figures  show  the  variation  in  calculated  boresight  error  as  a  function  of 
range.  These  results  reflect  a  wide  range  of  parameters  and  assumptions  about 
different  aspects  of  the  analysis.  First,  we  vary  km  and  the  PDF  of  the  surface 
heights.  For  both  these  cases,  the  results  shown  are  for  a  model  that  assumes 
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there  is  an  average  18.  5  dB  gain  over  the  3*  beamwidth  of  the  monopulse  differ¬ 
ence  pattern. 

In  Figure  9,  <Jg  is  plotted  versus  range  for  exponentially  distributed  surface 
heights  where  *  1. 5.  processing  losses  of  -S  dB  are  included,  and  the  effects 
of  shadowing  are  taken  Into  account.  Comparison  of  Figure  9  with  Figure  6  shows 
that  the  theoretical  boresight  errors,  Og  are  somewhat  less  than  the  experimental 
values. 


Figure  9.  Theoretical 
Calculations:  <3q  vs  Range, 
km  *  1.  S,  Shadowing, 
Exponential  PDF 


Figure  10  has  the  same  parameter  values  as  Figure  9.  except  that  k  ■  1.7, 

instead  of  km  ■>  1.  5.  Comparison  of  Figure  10  with  Figures  9  and  6  shows  that 

the  higher  value  for  km  results  in  slightly  poorer  agreement  between  theory  and 

experiment.  In  conclusion,  there  is  no  significant  effect  on  Cg  for  the  range  of 

k  values  used  here, 
m 


Figure  10,  Theoretical 
Calculations;  Jg  vs  Range, 
k_  »  1.7,  Shadowing 
Exponential  PDF 
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All  the  parameters  in  Figure  11  are  identical  with  those  in  Figure  9,  except 
that  the  PDF  for  the  surface  heights  is  Gaussian,  instead  of  exponential.  Compari¬ 
son  of  Figure  11  with  Figure  9  shows:  The  assumption  that  the  surface  heights 
are  all  normally  distributed  results  in  less  diffuse  multipath  power  entering  the 
receiver,  with  a  consequent  decrease  In  the  boresight  error. 

These  two  conditions  represent  the  extremes  of  the  possible  effects  of  surface 
height  distribution  on  boresight  error  for  these  two  PDF's.  The  initial  statistical 
analysis  of  the  terrain  heights  in  the  boxes  along  this  particular  trajectory  leads 
to  the  conclusion  that  the  exponential  PDF  was  most  appropriate  in  all  instances. 
Thus,  that  statistical  model  asserts  that  the  results  of  Figure  9  and  Figure  10  are 
the  ones  to  consider  rather  than  those  of  Figure  11.  Additional  statistical  approaches 
(including  the  multiple  observation  bivariate  formulation)  are  being  examined  to 
assess  their  conclusions  for  the  PDF's  of  these  same  regions. 


Figure  11.  Theoretical 
Calculations:  (Jq  vs  Range, 
k—  *•  1.5,  Shadowing. 
Normal  PDF 


The  next  effect  to  be  examined  is  the  inclusion  of  an  azimuthal  variation  In  the 
difference  channel  power  pattern  instead  of  an  average  value.  This  is  discussed 
in  Appendix  C.  Figure  12  shows  this  result.  Here,  km  ■  1.5,  the  surface  height 
PDF  is  assumed  to  be  exponential,  and  the  effects  of  shadowing  are  included.  Com¬ 
parison  of  Figure  9  with  Figure  12  shows  that  the  inclusion  of  the  azimuthal  varia¬ 
tion  In  power  pattern  results  in  a  decrease  in  diffuse  multipath  power  and  a  decrease 
in  boresight  error  Oq.  This  is  due  to  the  fact  that  the  width  of  the  glistening  sur¬ 
face  as  defined  in  Beckmann  and  Spizzichlno1  is  very  narrow  {approximately  40  m 
or  about  0.  1*  angular  extent).  The  calculation  thus  uses  an  average  difference 
pattern  gain  over  this  angular  region  •  20  dB  below  the  sum  channel  gain.  On  the 
other  hand,  in  the  approximation  that  there  is  no  azimuthal  variation  in  power 
pattern,  the  value  of  6  dB  below  the  sum  channel  gain  is  used  as  the  average.  This 
accounts  for  the  decrease  In  diffuse  power  and  boresight  error  when  the  actual 
azimuthal  variation  is  Included, 


Finally,  by  comparing  Figure  6  with  Figures  9  and  12.  we  conclude  that 
inclusion  of  the  azimuthal  variation  in  power  pattern  results  in  poorer  agreement 
between  theory  and  experiment  for  the  bores  ight  pointing  accuracy. 


S.4  Diffuse  Multipath  Fewer 

The  boresight  error  includes  a  diffuse  multipath  contribution  and,  for  complete¬ 
ness.  that  result  is  also  presented.  Figure  13  shows  the  variation  in  Pjjjpp  f°r 
the  same  parameter  values  used  in  the  case  of  Figure  12.  Comparison  of  Figure  13 
with  Figure  8  shows  that  the  diffuse  multipath  power  in  the  difference  channel  ia 
between  30  dB  to  50  dB  less  than  the  coherent  power  in  the  sum  channel. 


Figure  12.  Theoretical 
Calculations:  op  vs  Range, 
kjn  ■  1.5.  Shadowing. 
Exponential  PDF,  Azimuthal 
Pattern  Variation 


Figure  13.  Theoretical 
Calculations:  Diffuse  Power 

(PDIFF)  v*  Rang®'  km  *  *•  5* 
Shadowing,  Exponential  PDF, 
Azimuthal  Pattern  Variation 


5.5  Conclusion 

The  results  of  this  report  are  being  extended.  Improved  techniques  for  the 
estimation  of  statistical  parameters  characterizing  the  terrain  and  alternative 
methods  for  hypothesis  testing  of  the  PDF  of  the  terrain  heights  are  being  pursued. 
Further  improvements  involve  the  introduction  of  mean  surface  tilt  effects,  models 
that  contain  two  scales  of  surface  roughness,  and  the  azimuthal  variations  of  the 
normalized  rough  surface  cross  section,  0Q.  Preliminary  results  indicate  that 
inclusion  of  the  azimuthal  variation  of  a ^  results  in  excellent  agreement  between 
theory  and  experiment  for  boresight  pointing  error. 
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Appendix  A 

The  Shadowing  Function  for  Exponentially 
Distributed  Surface  H eights 


Brown10  has  derived  an  explicit  expression  for  the  shadowing  function  in  the 
case  of  backscatter  when  the  bivariate  PDF  for  the  surface  slopes  has  an  exponen¬ 
tial  form.  Here,  we  present  the  corresponding  results  for  bistatic  scattering, 
based  on  Brown's  work.  In  conformity  with  Eq.  (20)  in  Section  4.  1,  let  the  polar 
angles  of  incidence  and  scattering  be  denoted  by  9^  and  9g,  and  let  the  azimuthal 

angles  of  incidence  and  scattering  be  denoted  by  6.  and  0  ,  Also,  let  S  denote 
°  i  s  exp 

the  shadowing  function  when  the  bivariate  PDF  for  the  surface  slopes  is  exponentially 
distributed.  Then,  for  the  special  case  when  Qg  *  77  +  9^: 


S 


exp 


for  9  -  9  . 

s  i 


and 


S  =  r  1  ,  for  d  2  9. 
exp  +  1  s  l 


For  all  other  cases: 


exp  CQ  +  C,  +  1 
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For  normally  distributed  surface  heights,  the  expressions  for  CQ  and  C2  are  given 
by  Sancer. 9  Here,  for  exponentially  distributed  surface  slopes,  the  expressions 
for  CQ  and  C2  are  given  by 

Co  ‘  Ui/ff)  K2(xi>  *  1/2  +  <xi/2)  fKl(xi>Lo(lti)  +  Ll(xi)Ko(xi>! 
and 

C2  -  UJn)  K2(xa)  -  1/2  *  (xg/2>  (K1(xs)  Lq(x8)  +  L1(xa)  K^)] 
where 

x.  =  [yj6  / (2<7/T)J  cot  6  . 

\  =  tV®~ /<2a/X>J  cot  6g  , 

Kj(x)  =  Modified  Bessel  function  of  the  second  kind  of  order  i  , 

and 

Lj(x)  -  Modified  Struve  function  of  order  j  (see  Reference  20)  . 

20.  Abramowitz,  M. ,  and  Stegun,  1.  A.  (1988)  Handbook  of  Mathematical  Functions. 
U.S.  Govt.  Printing  Office. 
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Appendix  B 

Spsoulsr  Region  Characterization 

The  purpose  of  this  discussion  is  to  describe  the  region  contributing  to  the 
coherent  multipath  at  the  receiving  antenna.  For  the  case  of  interest,  the  angles 
of  incidence  are  relatively  small.  The  approximation  is  made  that  only  the  first 
Fresnel  zone  effect  is  significant;  successive  phase  differences  between  the  re¬ 
maining  zones  produce  terms  that  would  destructively  interfere  with  each  other 
at  the  receiver. 

The  Fresnel  zones  can  be  described  by  a  family  of  nested  ellipses.  The  curves 
In  Figure  3  show  the  range  dependence  of  the  zone  length,  width,  center  point,  and 
specular  point  location.  These  results  are  based  on  relations  found  in  Beckmann 
and  Splzzlchino.  * 

The  results  are  given  in  terms  of  target  distance  L  from  the  receiving  antenna 
along  the  ground  projection  of  the  target  trajectory.  Lj  is  the  distance  to  the 
specular  point,  h^  is  the  antenna  height  above  sea  level,  h^.  the  transmitter  height, 
and  ?8p  the  mean  height  of  the  subregion  containing  the  specular  point.  First, 


To  establish  the  Fresnel  zone  length  and  width,  we  define  the  path  length  distance 
relative  to  the  direct  ray,  as 
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‘"■a-S^t-V  (m, 


where 

6q,  hA>  and  1^  are  ail  much  smaller  than  L. 
Then,  the  center  of  the  elliptical  zone  XQ1  is  given  by 
~XL  ^  2(hA  ~  «,p> 


X01  *  (  t) 


X  L  ♦  (hA  +  -  2  z  )2 


the  semimajor  axis  intercept,  Xjj  is  given  by 


/x  l2\ 

(X  +  2  6  )1/2 

O 

XL+  <hA+  Kj.  -  2^/  _ 

(km)  , 


C<m), 


and  the  semiminor  axis  intercept,  V  j  j  is  given  by 


X  4*  2  0 

_ o 

XL  ,(hA  +  IVj.  -  2  Za/ 


T  1/2 


(m)  . 


It  should  be  noted  that  these  expressions  are  approximate  ones  and  that  there 
is  considerable  discussion  as  to  just  what  is  the  actual  region  that  contributes  to 
the  specular  multipath. 
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Appendix  C 

Tbs  Azimuth*!  Variation  of  tha  Monoputaa 
DHtantmPimm 


In  order  to  obtain  a  more  realistic  value  for  the  azimuthal  power  pattern  con¬ 
tribution  to  the  power  at  the  monopulse  receiving  antenna,  the  original  simple 
constant  value  was  replaced  by  values  based  on  the  actual  difference  (power)  pattern 
described  by  McGarty.  11  This  azimuthal  monopulse  difference  pattern  is  illus¬ 
trated  in  Figure  Cl. 

For  our  purpose,  the  intensity  (y-axis  of  the  figure)  was  converted  from  dB 
into  pure  numbers  (that  is,  -10  dB  =  0.  1,  -20  dB  =  0.  01,  and  so  on.  The  main  null 
in  the  difference  pattern  was  assumed  to  have  a  value  of  y  ■>  0  at  azimuth  angle, 
x  -  0.  0°  (x-axis  of  the  figure).  The  two  peaks  in  the  difference  pattern  that  occur 
at  x  =  ±2,5°  are  assumed  to  have  an  intensity  y  *  -4  dB  »  0.398.  Since  the  aver¬ 
age  sidelobe  level  is  less  than  or  equal  to  -45  dB  ■  0.0000318,  the  intensity  at 
x  -  5.  0°  was  set  equal  to  y  =  0. 

We  approximated  half  the  difference  pattern  (y  values  for  x  -  0)  by  fitting  a 
parabola  through  the  three  points  (x  =  0°,  y  =  0),  (x  =  2.  5®.  y  =  0,  398)  and 
(x  =  5®,  y  *  0.0).  This  parabola  Is  represented  by  the  equation: 

y  =  -0.  0637  x2  +  0.  3184x  . 

We  want  to  determine  the  relation  of  the  half  width  of  the  glistening  surface, 

£^s  as  defined  by  Beckmann  and  Spizzichino,  1  to  the  half  width  of  the  antenna 
pattern  footprint  on  the  ground  £pp.  Let  Rpp  be  the  range  from  the  monopulse 

antenna  to  the  point  on  the  ground  and  recall  that  xpp  =  5'  on  the  difference  pattern, 
the  value  beyond  which  the  intensity  y  -  0.  Then 

fpp  a  Rpp  tan  xpp  =  ^pp  *an  5®  , 

fGS  *  RFP  tan  XGS  * 
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In  the  computer  program  for  calculating  the  coherent  and  incoherent  power 
reaching  the  mono  pulse  receiving  antenna,  the  total  receiving  antenna  power  pattern 
is  taken  to  be  the  triple  product  of  the  acimuthal  pattern  given  by  Figure  Cl  (differ¬ 
ence  pattern)  multiplied  by  the  elevation  power  pattern  (difference  pattern,  given 
in  the  report  by  McCarty)1*  multiplied  by  the  sum  channel  gain  (82. 5  dB).  In  this 
Appendix,  we  are  considering  the  effects  only  of  the  acimuthal  power  pattern— the 
elevation  power  pattern  and  sum  channel  gain  have  been  normalised  to  unity  in  the 
present  form  of  the  expression.  Then,  the  average  power  PAV>  from  the  glisten¬ 
ing  surface  detected  In  the  difference  channel  is 


AV 


*GS 

J  y  (x)  dx 


or 


PAV  -  -0.02123  4s  +  0.159xgs 


{or  £gs  ~  * 


FP 


and 


P^y  *  -0.02 123  Xpp  0.  1S9  Xpp  for  f^g  2  fpp  • 

Note  that  in  the  last  expression  the  contributions  to  the  diffuse  power  from  pointB 
of  the  glistening  surface  beyond  fFp  are  considered  negligible  because  of  the  low 
antenna  side  lobes. 
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Figure  Cl.  Monopulse  Difference  Pattern 
in  Azimuthal  Plane 
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